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Abstract 

We analyze the asymptotic behavior corresponding to the arbitrary 
high conductivity of the heat in the thermoelectric devices. This work 
deals with a steady-state multidimensional thermistor problem, con- 
sidering the Joule effect and both spatial and temperature dependent 
transport coefficients under some real boundary conditions in accor- 
dance with the Seebeck-Peltier-Thomson cross-effects. Our first pur- 
pose is that the existence of a weak solution holds true under minimal 
assumptions on the data, as in particular nonsmooth domains. Two 
existence results are studied under different assumptions on the elec- 
trical conductivity. Their proofs are based on a fixed point argument, 
compactness methods, and existence and regularity theory for elliptic 
scalar equations. The second purpose is to show the existence of a 
limit model illustrating the asymptotic situation. 

Keywords: Seebeck-Peltier-Thomson effect; Joule effect; Nemytskii op- 
erator 
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1 Introduction 

In the recent years, the mathematical study of the thermistor problem has 
been of interest by various authors (see for instance [U EJ El EH [M] and the 
references therein) however the cross effects are neglected. Here we study a 
mathematical model for thermoelectric conductors, introduced in [10], tak- 
ing into account the presence of the Seebeck-Peltier-Thomson and the Joule 
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effects. Indeed, in the thermodynamics analysis, the Joule effect is given by 

|j| 2 /<,#) with ( C f. [23) 

j = -<r(.,0)(a(.,0)V0 + W), (1) 

representing the current density, 9 denoting the temperature, <fi is the electric 
potential, and the electrical conductivity a is a known positive function. The 
electrical conductivity is assumed temperature dependent and this is different 
at different places along the material due to the molecular structure. 

The Seebeck coefficient a is a given nonlinear function, dependent both in 
space and temperature, with constant sign observing that the sign of the See- 
beck coefficient corresponds to the sign of the Hall effect [22J. The Thomson 
effect is —da/dT(-,6)6'V6 • j, where d/dT means the derivative with respect 
to the real variable. Due to the first Kelvin relation, n(9) = 6ot{6), 9a(9) 
corresponds to the Peltier coefficient. Due to the second Kelvin relation, 
n{6) = 9a' (9), this coefficient is known as Thomson coefficient. Although \i 
is the only thermoelectric coefficient directly measurable for individual ma- 
terials [Hl[23], and the Seebeck coefficient appears as a{T) = J T ^-dt for 
some reference temperature T r , we keep the Seebeck coefficient as a given 
function as it is usual in the literature for the thermistor problem when this 
cross-effect coefficient is taken into account (2j El [33] . Thus, the resulting 
PDE's system is strongly coupled. 

The thermoelectric problem under study reads (see its derivation in [TC)| 



(V) Find the pair temperature-potential (9, 0) such that 

- V ■ {kV9) = a(-, 9)a(-, 9)(a(; 9) + — (-, 9)9)\ W| 2 + 

+<, 9)(2a(; 9) + ^(-, 9)9) V9 ■ V0 + <, 9)\ V0| 2 + g in O, (2) 

-V- (<r(-,0)V0) = V- (a(-,9)a(;9)V9) in ft, (3) 
W9-n=-a(-,9)9h on 9ft, (4) 
a(-,9)(V(j) + a(- 1 9)V9) -n = h on 9ft, (5) 

where ft is a convex bounded domain of IR n (n > 2), and it may represent 
electrically conductive rigid solids such as for instance thermistors, thermo- 
couples, resistive thermal devices (also called resistance temperature detec- 
tors) or thermoelectric coolers (see [3l [5j [35] and the references therein). 
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Notice that the boundary dQ is Lipschitz since every bounded convex open 
subset of M n always has a Lipschitz boundary [T5J Section 1.2]. 

Here n is the unit outward normal to the boundary dQ, g denotes the 
external heat sources and h denotes the surface current source. The thermal 
conductivity k is the known positive coefficient of the Fourier law. Finally, 
we remark that (jl]) is known as the linear Newton law of cooling with a(6)h 
representing the heat transfer coefficient. 

Our two main purposes are: (1) to improve the existence result stated 
in [10] for three dimensional space, i.e. the solution belongs to the Sobolev 
space W 1,P (Q) with p > n = 3. Indeed, we state the existence result for the 
multidimensional case since its proof is valid for any n > 2; (2) to pass to the 
limit on the thermal conductivity, k — > +oo, in order to show the existence 
of the limit model 



for some positive constant O solving an implicit scalar equation. The nonlo- 
cal problem ©-0 i s known as the shadow system, the heuristic designation 
introduced by Nishiura [31] in order to exhibit minimal dynamics displaying 
the mechanism of basic pattern formation. Here the shadow system turns 
out the dynamic relation among the trivial rest states (constant solutions) 
and the large amplitude voltages. The present result generalizes the homoge- 
neous Neumann boundary value problem already studied for nonlocal elliptic 
problems [7J. 

The regularity assumptions about the domain could be weakened if dif- 
ferent techniques are provided (see [U [13], [HJ [T7J [281 ES] and the references 
therein). We refer to [12] the existence of singularities of electromagnetic 
fields at corners and edges of a bounded Lipschitz domain with piecewise 
plane boundary. 

The contents of this work follow. In next Section we state the as- 

sumptions and the main results. Section [3] deals with existence and regularity 
results for auxiliary problems. We prove in Section H] the existence result for 
n > 2 when the electrical conductivity is assumed be a uniformly continuous 
function, and we obtain in Section its limit model. The proof of the exis- 
tence result valid for n = 2 when the electrical conductivity is assumed be 
discontinuous on the space variable and the corresponding asymptotic limit 
model are postponed in Section [61 



V- (a(-,6)V0) = infi; 
a(-,Q)V(/>-n = h onffi, 



(6) 
(7) 
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2 Assumptions and main results 

In order to establish the existence results we assume the following set of 
hypotheses on the data. 

(HI) k is a positive constant. 

(H2) (j:f2xIR— )> IR is a Caratheodory function, i.e. measurable with respect 
to x G Vl and continuous with respect to T G M, and furthermore 

3o#, o* > : a # < T) < a*, VT G R, a.e. in Q. (8) 

(H3) <? belongs to the Lebesgue space L P//2 (f2), p > 2, and h G C(dQ) is 
such that verifies the compatibility condition 

/ h ds = 0, (9) 

where ds represents the element of surface area. 

(H4) a G C{Q x R) is such that 

3L a > : |a(xi,Ti) -a(x 2 ,T 2 )| < L a (\x l —x 2 \ + |Ti -T 2 |), (10) 

for all Xi,X2 G and for all T 1? T 2 G R, with |x| representing the 
euclidean norm and |T| the absolute value of a real number. Moreover, 
for all x G f2, the mapping T h-> a(x,T) is continuously differentiable 
in M and its derivative satisfies 



3,*>0: |_ (I ,D| < | ^ |n if j T j; L (11) 

The following two different cases will be addressed: 
(H4) + for materials with positive a: 

3a#, a* > : a # < a(x, T) < a*, \/x G Q, VT G R. (12) 

Moreover, there exist an open subset r C dQ such that meas(r) > 
and meas(<9f2 \ f ) > and constants h# > and hi < such 
that 

mm^M Mx) >(J#. a-e-^r 
Cia # I Aii, a.e. a; G S, 



where S := dQ \ T and C\ denotes the continuity constant of the 
embedding H\Q) e — >■ L 2 (£), i.e. 

W, E <Ci(l|VCn+l|f|llr) (14) 

for every 9 G H 1 ^); 
(H4)_ for materials with negative a: 

3a # , a* > : -a* < a(x, T) < -a # , Vx G Q, VT G E. 

(15) 

Moreover, there exist an open subset T C <9fi such that meas(r) > 
and meas(<9fi \ f) > and constants h# > and foi < such 
that 

> _ min{fc,«jft # } < j -J #J a.e. * G T 

da* [ -/ii, a.e. x G L. 

Remark 2.1 In particular, [Lfy holds for the unity function then we have 

|£|<Ci|r|, (17) 

where \T\ =meas(T) and |S| =meas(dQ \ T). 
Some remarks on the assumptions. 

1. The heat conductivity can be an uniformly continuous function on both 
variables verifying 

3k # , k* > : k# < k(x, T) < k#, \/x e SI, VT G K. 

Since our purpose is to study the asymptotic behavior as k — > oo, we 
assume it as constant. 

2. The case T = dfl is excluded, from the fact that the Gauss theorem 
yields the necessary condition ([9]) of the existence of a solenoidal func- 
tion satisfying (P) and © (cf. ©). 

3. From the assumption flllD we obtain that |f^(x, T)\ < for all x G Q 
and for all T 6 K, and that the Thomson coefficient is bounded, i.e. 
\fi(x,T) \ = \T^(x,T)\ < n#, for all x G Q. and for all T G R. 
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4. By the Weierstrass Theorem, any continuous function defined on a 
compact set (of M n+1 ) is bounded. Then the upper bound in (IT21 as 
well as the lower bound in f|T5|) could be given reduced to \T\ > 1. 

We define the Nemytskii operators 

(7(6) = a(;9(.)); a(9) = a(;9(.)); ^(9) = ^M(-)), 

that map L l (Q) into L q (Q), for all q < oo. Their designation is kept in order 
to clarify the presentation. 

We endow the Sobolev space W 1,P (Q), p > 1, with the equivalent norms: 

• for the temperature solution 

||0||i,p,n = ||V0|| P) fi + ||^||p,r; 

• for the potential solution 

U\k P ,a = HV0IU, 
considering the correspondent Poincare inequality. 

Definition 2.1 We say that (9, <p) is a weak solution to (E))-([3J) if (9, <p) e 

W 1,P (Q) 2 , for p > n, and it satisfies 



k / V9-\7r]dx+ / a(9)h9r]ds 
Jn Jan 



o{9) ( a(9)(a(9) + ^{9)9)\V9\ 2 + (2a(9) + ^{9)9)V9 • V^dx 



+ [ (a(9)\V(p\ 2 + g)r]dx, V77 £ (18) 

a{9)V(p ■ Vrjdx = 



= - [ a{9)a{9)V9 ■Vr 1 dx + [ hr]ds, V77 G VK 1,p '(f2), (19) 
where p' denotes the conjugate exponent to p, p' — p/(p — 1). 
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Remark 2.2 If Q, is of class C 0,1 and mp > n, then the Morrey-Sobolev 
embedding holds 

W m ' p (Q) c°' m - n/p (n). 

The vector field n belongs only to L°°(<9f2) if it is the unit outward nor- 
mal vector to the boundary of C 0,1 domains. When p > n, the embedding 
W l ' p '{yt) m> Lv n /tP n - n -P\VL) is valid. Thus the L p / 2 behavior of the quadratic 
terms |V6*| 2 7 V6- V</>, |V0| 2 is meaningful on the right hand side of (TJ|) since 
p/2 > pn/(p + n). 

Let us extend the existence results whose can be found in [10]. The first 
main theorem states the existence of weak solutions to the problem under 
study, strengthening the assumption (H2), i.e. strengthening the regularity 
on o. 

Theorem 2.1 Assume n > 2, (HI) and (H3)-(H4) hold. Additionally we 
assume that a : x IR — > M, is a uniformly continuous function satisfying ([#)) 
and the smallness condition |53|) is satisfied. Then the variational problem 
(Q2P-([Zg; admits a weak solution (6,<j>) G W 2 ' pn ^ p+n \n) x W X *(P), with 
p > n, in the sense of Definition HT7| such that J dn <j)ds. Moreover, the 
following estimates hold 

||V0|| 2 , n < — (o-*a*\\Ve\\ 2 ^ + C 2 \\h\\ 2tdn ); (20) 

||0||i,p,n < C p ; (21) 
fc||Ve||| n < C 3 (a # M<^ + 

+(AC 2 + BC P \\ V0|| 2 , n + a*\\ V0||^ + \\g\\ p/ 2,a)C p ) ,(22) 

where C 2 denotes the Poincare-Sobolev continuity constant of the embed- 
ding H 1 ^) 7- L 2 (dQ), C3 denotes the continuity constant of the embedding 

w^isi) c(n), 

A: =( T # a # (a # + /i # ), 5 := a # (2a # + /i # ); (23) 
and C p denotes a positive constant independent on k if k > a#h#. 

If we assume a G C(Q x R) as in [10] . the Nemytskii operator maps C(Q) 
into C(f2) which implies the uniform continuity on the spatial variable. 

Next, we establish the existence of a solution (f> = 0(0) to ©-([Tj) where 
G is solution of an implicit scalar equation. 
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Theorem 2.2 Under the assumptions of Theorem \2.1\ there exist weak so- 
lutions (0 fc ,0fc) G W 2 'P n /( p+n )(fi) x W l > p {vi), withp > n, to the variational 
problem / f7^] )-/ f7Pj) . such that 

9 k ^Q in H\tt), <j) k ^(f) in H 1 ^), 

with 6 solving the implicit scalar equation 

6 I a(-,Q)hds= [ o-(-,Q)\V4>\ 2 dx+ [ gdx (24) 
Jan Jn Jn 

and (j) solving (0|)-([?|). Moreover 

°#L\ V( t ) \ 2dx + !n9 dx < @ < Jn\^ ( t ) \ 2dx + Jn9 dx ^ 
a* J dQ \h\ds ~ ~ a#h#\T\ + a#hi\T,\ 

REMARK 2.3 In face of if the Seebeck coefficient is only a function on 
the temperature then (24\ ) reads 

0= / a(-,B)\V(p\ 2 dx+ [ gdx. 
Jn Jn 

The presence of a generic heat source invalids the limit model |6f-([2|) con- 
nected with the original thermoelectric problem introduced in FTO^ . 

The above results can be proved if the convexity of Q is replaced by weaker 
assumptions, for instance when Q is a plane bounded domain with Lipschitz 
and piecewise C 2 boundary whose angles are all convex [T9l p. 151], or when 
Q is a plane bounded domain with curvilinear polygonal C 1,1 boundary whose 
angles are all strictly convex [T91 p. 174]. For general bounded domains Q of 
M n with Lipschitz boundary, it is known that the integrability exponents for 
the gradients of the potential and temperature solutions may be larger than 
3 [T71 [36] , if the restriction to the case of uniformly continuous coefficients in 
(I18p - (ll9p is assured. However, a generalization for such nonsmooth domains 
of Theorem 12.11 and its limit model is not a direct consequence. Indeed, new 
proofs will be needed because the compact embedding W 2,pn ^ p+n \Vt) ' > r > 
W l,p (Q) is crucial to provide the weak continuity of the operator in the fixed 
point argument. 

For the two-dimensional limit model, let us show the existence result 
under the minimal regularity on o. 
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Theorem 2.3 Assume n = 2 and (H1)-(H4) hold. Then the problem (0^-([3P 
has a weak solution in the sense of Definition ^. 1\ for some p > 2 = n, under 
sufficiently small data. 

Proposition 2.1 Under the assumptions of Theorem \2.3L there exists <ft solv- 
ing 1?]) with solving the implicit scalar equation ( f^p . 

The study of the existence of three-dimensional weak solutions to the 
variational problem ( |T5|) - (fT^|) . under the assumption that the mapping x G 
Q i — y o~(x,T) is discontinuous for every T £ R, is still an open problem. 
When the coefficient o of the principal part of the divergence form elliptic 
equation (J3]) is a discontinuous function on the spatial variable, it invalidates 
the smoothness of the solution as is carried out in the literature [6j |2U l26| 

1231132]. 

3 Auxiliary results 

The existence of a unique solution G such that J da (f)ds = 0, to an 

auxiliary problem is consequence of Lax-Milgram lemma, for details see 

Proposition 3.1 ([10, Theorem 4.1]) Let the assumptions (H2) and ( fjjj) 
or ([73]) be fulfilled. Assume that n > 2, Q G C ' 1 , £ G H\tt) and h G L p (dtt) 
verify (G|) for p > 2(n — l)/n. Then there exists a unique weak solution 
(f) G H 1 ^), such that J dQ (j)ds = 0, to the variational problem 

/ <r(£)V0 • Vrjdx = - <7(£K0V£ • Vr?cfe + / ^ds, (26) 
Jn Jn Jan 

for all rj G H 1 ^) and in particular for all rj G W 1,p ' (Vl). 

Next, we establish some regularity for the potential auxiliary solution. 

Proposition 3.2 Let p > n, h G L p (dVt) verify (TjJ), £ G W 2 ' 9 ^) with 
q = pn/(p + n), (H4) be fulfilled, and <p G H 1 ^) solve the problem k2h^) . 
If o : Q x R — )■ R zs a uniformly continuous function verifying (0|) ; i/ien 
V0 G L p (fi) and verifies 

||V0|| Pj n < C7(n,p,n,a #) (7#)(a#a#||VeiUn+ ll/ilUm). (27) 
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Proof. For p > n, we have q > n/2, 

£ G W 2 ' q (Q) m- C°' 2 - n/9 (fi), a(OV^ G W ll9 (fi) m> L g(n - 1)/(n - g) (9fi), 
a(f)Vf • n G L g(n - 1)/{n - 9) (9fi) = L p(n - 1)/n (<9ft) iy~ 1/p ' p (dfi), 

for n G L°°(<9ft). Moreover, /i/<r(0 G W^^an) and V£ G L P (Q) implies 
that V- (ex (£)«(£) V£) G (W 1,p '(f2)) forp>n. By appealing to the regularity 
theory [TSJ, [TH] for the solution <fi G i/ 1 (Q) of the boundary value problem (in 
the sense of distributions) 

V-(<r(OV0 + <r(Oa(OVO = O in fi 
^(O(V0 + a(OVO-n = /i onffi, 

and observing that £ G W 1,p (fi) ^ C 0,1 ~ n / p (f2) for p > n warrants that 
cr(-,£) is uniformly continuous, then the regularity of weak solutions relative 
to W hp (Q) and the estimate (F27D arise. □ 

The following result deals with the existence and uniqueness of a strong 
temperature auxiliary solution. 

Proposition 3.3 Let p > n, £ G W^iVt), (HI) and (H3)-(H4) be fulfilled, 
and <f) G H 1 ^) solve the problem If2h]) . If a : Q x R — )■ R zs a uniformly 
continuous function verifying then there exists a unique weak solution 
e W 2,im/(p+n)(ty solving the problem, for all r] G 

fc / V6-Vr)dx + / a{C)hdr 1 ds= [ (F(-, £, V£, V0) + #) r/dx, (28) 
Jan Jn 

with F : Q x R 2n+1 -»■ R de/med as F(x, T, a, b) = 

a(x, T) (a(x, T)(a(x, T)+^(x, T)T)\a\ 2 +(2a(x, T)+^(x, T)T)a-b+|b| 2 ) . 

Proof. The existence and uniqueness of G W 1,P (Q) is consequence of 
Propositions 13.11 and 13.21 By appealing to the elliptic equations theory [19], 
from F(£, V£, V</>) + g G L p / 2 (f2), the regularity theory for the Laplace 
equation in convex domains guarantees the existence of a unique solution 
6 G W 2 ' pn ^ p+n \n) of the Robin problem 

-kA9 = F(-,Z,VZ,V0) + g in ft; 
kV6 • n + a(£)M = on dfi, 
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taking into account that the Korn perturbation method rT9J pp. 107-109] can 
be adapted if the coefficient a(-,£)h G C{dVL) is such that the assumption 
f fTOj) holds. For this, we observe that £ G iy 1,p (fi) C(fl) and we recall 
(H3)-(H4). □ 

For the regularity of the potential auxiliary solution when it is the 
unique weak solution for Neumann problem to an elliptic second order equa- 
tion in divergence form with bounded and measurable coefficient, we can 
prove the following result. 

Proposition 3.4 If the assumptions of Proposition \3J\ are fulfilled with p = 
2 and Q is convex, then there exists a constant e > such that the weak 
solution G F^fi) of (ESP belongs to W l,2+e (U), i.e. 

||V0|| 2+e ,n < K 2 {a*a*\\Vi\\ 2 ,n + \\h\\ 2 ,on), (29) 
with a constant K 2 > only dependent on the data. 
Proof. Denote the operator A by 

(Acf), 7]) = / <r(0Vf Vrjdx. 
Jn 

Then G H 1 ^) is a weak solution to the second order elliptic differential 
equation Au = F, under 

F = V • (cr(0a(0V0 + h G (tf 1 ^))' ^ (W^(fi))', Vp > 2. 
Since the boundedness property 

is fulfilled, considering that £ G and the assumption (jSj) on a holds, 

then the Neumann version of the general result on the higher regularity for 
weak solutions to the mixed boundary value problems (cf. [20l Theorem 1], 
also [21]) guarantees that G W 1,2+e (£i) for some e > 0. □ 

Although Proposition 13.41 is valid for any dimensional space (n > 2), 
we only used it for n = 2. Let us precise its application in the following 
proposition. 

Proposition 3.5 Let £ G W 1,2+e (fi), (H1)-(H4) be fulfilled, and G H\Q) 
solve the problem (fgfip. If the assumptions of Proposition^^ hold, then there 



exists a unique weak solution 9 G W 2,2p ^ p+2 \Vt) solving (TJ2|) with p = 2 + e. 
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Proof. The imperative requirement of the embedding W 1,2+e (Q) <^-> C(Q) 
yields that a(-,^)h G C{dVL) provided by (H3)-(H4). Thus, Proposition 
13.41 ensures that the argument of the proof of Proposition 13.31 is still valid, 
concluding the claim. □ 

4 Proof of Theorem 12.11 

First we recall the Tychonoff extension to weak topologies of the Schauder 
fixed point theorem [TBI PP- 453-456 and 470]. 

Theorem 4.1 Let K be a nonempty compact convex subset of a locally con- 
vex space X . Let £ : K — )■ K be a continuous operator. Then C has at least 
one fixed point. 

If we provide any Banach space with the weak topology, every closed ball 
is convex and weakly sequential compact. 

In order to apply Theorem I4.1[ let us consider the operator C defined in 
a closed ball B R C W 2 ' pn ^ p+n \n) such that 

£:(G5 fi 4^0G W/ 2 ' pn/(p+n) (^), 

where G W 1,p {Vt) solves the problem ( 12"B"|) . for all rj G W 1,p (£1), and 9 solves 
the problem (Kg)) . 

The existence of a unique solution G W 1,P (Q), such that f m <frds = 0, 
to the problem ( )26|) is consequence of Propositions 13 . 1 1 and 13 . 2 [ and it verifies 

(EZD. 

Hence, for p > n, we find 9 G W 1,P (Q) from Proposition 13.3} and the 
estimate 

||^|| 2 W(P+"),n < ^(^llVeilJ^ + ^llVeiUollV^II^ + ^llV^II^+lbllp/^n), 

(30) 

is verified with K denoting a constant dependent on Q, n and p, A and B 
given by (123"]). and x := min{A;, a#h#} + a*hiC\ > 0. 
Thus C is well defined. 

Next, let us prove that C(Br) C -Br. Let £ G -Br be arbitrary and (0, 0) 
be the corresponding solution solving ( 126|) and ( 128]) . Thus (127)1 and (!30|) read 

||0||i, P ,n < C{a*a*R+ \\ h\\ p , m ); (31) 
\\0h P n/( P+ n),n < K(AR 2 + BR\\V<t>\\ p ,n + V0||^ + |M| pA n), (32) 
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with C = C(n,p,Q,a#,a # ). Inserting ( 13"Tj) into ( 13"2"j) it follows 

||#||2,pn/(p+n),c < a 2 R 2 + aiR + a , 

where 

a 2 = Ka#a*(l + Ca#)(a*(2 + Ca*)+ I 2*); 
ai = ACV # (2a # (l + a # )+/i # ) \\h\\ PtdQ ; 
a = A" (CV # ||/i||p )an + \\g\\ P /2,n) ■ 

Therefore, £(£) = G -Br if and only if a 2 R 2 + (aj — + cto < 0, i.e. for 
instance if the smallness condition 

ai < 1 and 4a a2 < (1 — ai) 2 (33) 

is assumed. 

4.1 The weak sequential continuity of C 

Let {( m } m6 R be a sequence in B R verifying 

£ m in W 2 > pn/(p+n \tt) ^ x ' p (fi), (34) 

and (</> m , # m ) is the correspondent solution to fl26|) and (1281 . for each m e N. 
From the estimates fl27|) and fl30l) we can extract a subsequence, still labeled 

by (<f) m ,0m), such that 

m ->> in W^fi), # m ->> in W 2 ' pn/{p+n \Q) W lj, (Q). 

Thanks to Remark 12.21 it follows 

U^i, m ^9, m ^0 in C°' l - n l p {n). (35) 

In particular, J dQ (j) m ds = 0—7- j gQ <pds = 0. By the continuity of the Ne- 
mytskii operators a and a, we can pass to the limit in as m tends to 

infinity, concluding that <fr G W l,p (Q) is the limit solution, i.e. it verifies (|26p . 

In the sequel, let us pass to the limit in ([28]) m as m tends to infinity. 
First, the mapping £ e L 1 (fi) h-> a(£) G £ r (fi), for all r < +oo, is contin- 
uous by (H4), thus the passage to the limit to the left hand side of (|28|) is 



13 



straightforward. In order to study the RHS, we define 



I 



I. 



2,m 



:rn 



L 
J 

J 

Jn 



0-(Cm)a(Cm)(a;(Cm) + ^(£m)£m)|V£ m |Vfe; 

dot 

0-(Cm)(2a(Cm) + 7^(£m)£m)V '£ m • V 4> m Vdx] 
c{{,m)\V4> m \ 2 VdX- 



Recalling Remark Owe get rj G W 1 * '(fi) ^ xjm/dm-n-p) ^ ^ LP/(?>- 2 )(fi) 
for p > n. 

From (I54"|) we have |V£ m | 2 — )■ |V£| 2 in L P//2 (f2). Considering that the 
mapping £ G i-)- er(£)a; 2 (£)r7 e L P// ( p_2 ' is continuous thus the 

first term in passes to the limit as m tends to infinity. Using (134"]) and 
(135]) we have £ m |V£ m | 2 — > £|V£| 2 in L P//2 (f2). Considering that the mapping 
£ G L\tt) ^ a(0«(0|f (0»7 G £ p/(p ~ 2) ( fi ) is continuous thus the second 
term in J 1>m passes to the limit as m tends to infinity. 

Analogously, we take to the limit in J 2jm observing that the strong-weak 
convergence product V£ m • V0 m — ^ V£ • V0 holds in W^iVt). 

In order to be in conditions for finding that 9 is a limit solution, let us 
prove the continuity of the solution mapping £ G W 1,P (Q) h-> = 0(£) G 
iy 1,s (f2) in the strong topology for s = 2pn/ (p + n) < p. Take the difference 
of (1261) ™ and ( 1261) verified by the solutions m and 0, respectively, and choose 
V — 4>m — (ft as a test function. Thus, it results 



Then, we conclude that V0 m — )■ V0 in L 2 (f2), and consequently V0 m —> V0 



a.e. in ft and |V0 m | 2 -> |V0| 2 in L s / 2 {Vt) ^ L pn /( p+n )(fi). Thus J 3>m passes 



to the limit as m tends to infinity, concluding the proof of weak continuity 
of the operator C 

Then the Schauder fixed point theorem can be used and it guarantees the 
existence of (9, 0) in the conditions to Theorem 12.11 



C7#||V(0 m - 0)||| n < / ((7(0 " ^(£m))V0 • V(0 m - 0)cfe + 
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4.2 The validation of the estimates fl2~QT)-fl2~2~T) 

Let (9,(f>) G W 2 > pn/i - p + n) (Q) x W^(n) be a weak solution to the variational 
problem ffT8|) -ffT9|). 

Choose rj = (ft G W 1,P (Q) as a test function in f|T9|) . Using (JS}, the upper 
bound of | a | and the Sobolev-Poincare inequality then f[2"Uj) holds. 

From the regularity theory for the Robin-Laplace problem and by virtue 
of the existence of a solution 9 G W 2, ' pn ^ p+n \Vl) we proceed as in ( 130|) now 
for > ot#h# resulting the estimate 

||W|U < K{A\\V6\\l n + S||V0|UI|V0|Un + a*\\V<t>\\l n + IMUn), 

with K denoting a constant independent on k. Combining this result with 
the estimate (1271) with £ replaced by # and using (1331) we conclude (12T|) . 

Choose 77 = # G W /1,p (fi) as a test function in (1181) . Then applying the 
Holder inequality and using the assumptions (H1)-(H4) it follows 

k\\V8\\ln < « # M^2,£ + (^l|V0||I jn + £||W|| 2 , n ||V0|| 2 ,n + 

+^ # l|V0||^ + |M|i,n)||0||oo,n- 

This yields the estimate (122|) . 



5 Proof of Theorem 12.2 

For each given k > 0, let (9k, 4>k) be a solution to (fl8l) - (fr9j) in accordance 
with Theorem 12.11 From estimates (|22|) and (|20|) there exist subsequences 
still denoted by 6^ and <pk such that, for k — » +oo, 

V0 fc -> in L 2 (fi); ^ ^ 6 in W 1,p (fi) M-M- C(O); 

fc - in H\Q), 

with O constant on Q. Hence we can pass to the limit in ffT9]) as k tends to 
infinity resulting 

[ a(;e)V<j>-Vr)dx= I hr]ds, V77 G W 1,p ' (£1), 
Jq Jan 

or equivalently ©-([7]). 
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In particular, if we take 77 = 1 in (1181) we can pass to the limit as k tends 
to infinity resulting ([23]). 
Using (JS} it follows 

a # / \V(j)\ 2 dx+ gdx<Q a(-,Q)hds < a* / \V(j)\ 2 dx + I gdx. 
Jn Jn Jan Jn Jn 

Taking into account that the assumption (H4) + or (H4)_ and also ( ITT]) imply 

/ a(-,Q)hds < a* \h\ds; 
Jan Jan 

I a(-,e)hds > a # h # \T\ > 0, 

Jan 

we derive ([2I3]) , concluding the proof of Theorem 12.21 

6 The two-dimensional limit model 
6.1 Proof of Theorem Q 



Arguing as in Theorem 12.1] the Schauder fixed point argument can be ap- 
plied. For this, it is sufficient to see that the regularity relative to = 
W 1,P (Q) (Q bounded) can be applied for the unique weak solution of the vari- 
ational problem (126]) for p = 2 + e > 2 = n in accordance to Proposition 13.41 
Thus Proposition 13.51 guarantees the existence of 9 £ W 2 ' 2pj/(j>+2 \il) verifying 
(130]) . For every f e B R , inserting (J29]) into (]30]) it follows 

\\eh,2 P /(p+2),n < K(AR 2 + BRK 2 {a*a*R+\\h\\ 2m ) + 
+a*K 2 2 (a*a*R + \\h\\ 2 , dn ) 2 + \\g\\ p/2 ,n). 

Next arguing as in Section H] it leads to a smallness condition. From the 
continuity of the Nemytskii operator a due to the Krasnoselski Theorem we 
can proceed as in Section 14.1] considering that I^m passes to the limit as m 
tends to infinity since 

°{U) ^ weakly* in 

|V0 m | 2 r/ -)■ |V0| 2 r7 inL 1 ^). 

This concludes that we are in the conditions of applying Theorem 14.11 
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6.2 Proof of Proposition 12.11 



We can proceed as in Section [5] considering the existence of the sequence of 
solutions is provided by Theorem 12.31 and the proposition follows. 
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